Fourier-transform spectroscopy of Sr 2 and revised ground state potential 
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Precise potentials for the ground state X 1 S^" and the minimum region of the excited state 2 E J of 
Sr2 are derived by high resolution Fourier-transform spectroscopy of fluorescence progressions from 
single frequency laser excitation of Sr2 produced in a heat pipe at 950 °C. A change of the rotational 
assignment by four units compared to an earlier work (G. Gerber, R. Moller, and H. Schneider, J. 
Chem. Phys. 81, 1538 (1984)) is needed for a consistent description leading to a significant shift 
of the potentials towards longer inter atomic distances. The huge amount of ground state data 
derived for the three different isotopomers Sr2, Sr Sr and Sr Sr (almost 60% of all excisting 
bound rovibrational ground state levels for the isotopomer 88 Sr2) fixes this assignment undoubtedly. 
The presented ground state potential is derived from the observed transitions for the radial region 
from 4 to 11 A (9 cm -1 below the asymptote) and is extended to the longe range region by the 
use of theoretical dispersion coefficients together with already available photoassociation data. New 
estimations of the scattering lengths for the complete set of isotopic combinations are derived by 
mass scaling with the derived potential. The data set for the excited state 2 1 E^ was sufficient to 
derive a potential energy curve around the minimum. 

PACS numbers: 31.50.Bc, 33.20.Kf, 33.20.Vq, 33.50.Dq 
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I. INTRODUCTION 

Over the past decades the Sr2 molecule has been the 
subject of many different studies. In 1977 Miller [2,2 dis- 
covered by the analysis of matrix isolation spectroscopy 
two excited 1 E+ states belonging to the (5 1 S 4- 4 1 D) 
and the (5 1 S + 5 1 ?) asymptotes and investigated the 
ground state X 1 S+ A fewyears later, in 1980 and 1984, 
Bergeman [3J and Gerber [4j performed heat pipe exper- 
iments and studied the transition from the ground state 
to the 2 1 £+ state (5 : S + 5*P asymptote), which they 
called A state. In 1992 Bordas @ discovered by deple- 
tion spectroscopy on a molecular beam a higher lying 
1 Ti u state, which was assigned in 1996 by ab initio calcu- 
lations [6j as the 3 1 H U state correlating to the asymptote 
5 3 P+5 3 P. Several other ab initio calculations were done 
on Sr 2 @, Is H HO, Qj| ■ In the most recent calculation 
[ll| the basis set has been chosen to reproduce best the 
ground state potential reported by Gerber J4|. Since for 
Strontium the (5 1 S + 4 X D) asymptote lies below the (5 1 S 
+ 5^) and thus the A 1 E+ state is not the lowest X I]+ 
state, we renamed this state to 2 1 E+ to avoid any fur- 
ther confusion in this work, according to the ab initio 
calculation @. 

Currently there is high interest on ultracold ensembles 
of Strontium atoms and high precision spectroscopy on 
Strontium, because it could be a candidate for an opti- 
cal frequency standard [13 . fl3l UA . fl5j . Very recently, 
Zelevinsky [li| proposed to measure precisely the time 
variation of the electron-proton mass ratio by the use 
of ultracold Si'2 molecules trapped inside an optical lat- 
tice. A lot of trap experiments [lj, [l8|, [l9|, [20] at cold 
and ultracold temperatures are reported for which the 
reliable knowledge of the interaction properties like scat- 
tering lengths of the different isotopes in different elec- 
tronic states would be of advantage. First estimates of 



the Strontium scattering lengths for homonuclear colli- 
sions of the isotopes 88 Sr and 86 Sr in their electronic 
ground states are already derived from photoassociation 
data [2l[. But there is currently no spectroscopic work 
available, which is sufficiently precise to derive such quan- 
tities from the scattering wave-function by directly solv- 
ing the radial Schroedinger-equation with a complete in- 
teraction potential. This would also yield the precise 
scattering lengths of different isotopes and isotopic com- 
positions, if the Born-Oppenheimer approximation holds 
for this molecule. 

Our current goal is to derive such precise potentials 
including the long-range part. We started our work by 
using the same excitation path as in the earlier measure- 
ments (4[ on the 2 1 E+ -- X X S+ transition, where we 
recorded 37 very dense fluorescence spectra with many 
overlapping bands. For their analysis we applied a newly 
written automated software to accelerate the identifica- 
tion and to reduce the probability of transfer error in 
producing data files for the later fitting routines. 

The paper is divided into the following sections: Sec- 
tion [TT] describes the experimental methods and the ap- 
paratus, section ITTT1 gives the models for data reduction, 
section ITVl describes how the analysis and the assignment 
was done and how the automated software works, section 
IVl presents the resulting potentials and Dunham coeffi- 
cient sets and describes how they were produced and fi- 
nally section [VI] discusses the usability of the presented 
data including a short outlook. 



II. EXPERIMENT 

The applied experimental methods are similar to the 
ones in [221 ] for Ca2. To form a gas sample of Sr2 
molecules, about 10 g of Strontium are heated inside a 



stainless steel heat pipe to a temperature of 950 °C. As 
buffer gas 20 mbar of Argon is introduced. The windows 
at the ends of the heat pipe are on both sides broad band 
antireflection coated. 

A special difficulty of the work with alkaline earth met- 
als like Strontium is the relatively high atomic vapor pres- 
sure of a few mbar at the melting point of 777 °C, which 
prevents good heat pipe operation conditions. Thus a 
significant part of the Strontium vapor condenses in re- 
gions of the heat pipe where the temperature is below 
the melting point and crystals grow towards the center 
and so towards the path of the laser beam. The time 
these crystals need to reach the laser beam path strongly 
depends on the temperature and the buffer gas pres- 
sure. Higher temperatures, which would lead to higher 
production rates of molecules and thereby stronger sig- 
nals, would also severely shorten the possible measure- 
ment time. This can be partially compensated by higher 
buffer gas pressures, but an increase of the pressure does 
not only rise the magnitude of pressure broadening and 
shift, but also leads to additional new problems. Every 
time Strontium is melting somewhere in the heat pipe 
or the buffer gas pressure is changed, some kind of fog 
occurs near the ends of the heat pipe, which mainly con- 
sists of fine drops of liquid Strontium, as it became clear 
when the fog accidentally reached a window and instantly 
turned it into a metallic layer. This fog produces strong 
laser stray light, increasing the noise level in the spectra, 
until it disappears by waiting. The waiting time and the 
fog density strongly increases with the buffer gas pres- 
sure. 

At a temperature of 950 °C and 20 mbar buffer gas it 
is possible to measure for two to three hours until the 
Strontium crystals block the laser beam path. When 
that occurred the heat pipe was shifted inside the oven 
to heat the ends and melt away the Strontium crystals. 
By this procedure most of the liquid Strontium gets ad- 
sorbed to a mesh and is transported back to the center 
of the heat pipe, but a smaller part is also transported 
further to the ends and reaches the windows to a small 
percentage. Though we did not achieve proper conditions 
of heat pipe operation where permanent recycling of the 
material would take place, it was nevertheless possible to 
record all spectra analyzed for this work without refilling 
and cleaning the heat pipe. 

The laser excitation is done by a Coherent CR 699 ring 
dye laser operated with rhodamine 6G in the frequency 
range from 17000 to 17600 cm -1 . The fluorescence light 
emitted antiparallel to the direction of the laser beam 
is imaged into a Fourier-transform spectrometer (Bruker 
IFS 120 HR), which is used with a resolution of 0.05 
cm -1 . Most of the recorded spectra are averaged over 10 
scans, in few cases also up to 100 scans are taken. 

Figure [1] shows a typical spectrum. As one can see the 
density of lines and the amount of excited progressions 
under single frequency laser excitation (width 20 MHz) 
is relatively high, resulting from the large reduced mass 
of the Strontium dimer and the high temperature, which 



leads to a significant thermal population of all existing 
rovibrational levels of the flat ground state potential of 
an excimer type molecule like Sr2. 



III. POTENTIAL ENERGY CURVES AND 
DUNHAM COEFFICIENTS 

The rovibrational energies of electronic states of 
type 1 S are calculated by the simplest one-dimensional 
Schrodinger equation, which includes the radial kinetic 
energy, the potential energy curve (PEC) V(R), and the 
centrifugal energy for the rotational state J. V(R) is of- 
ten called the Born-Oppenheimer potential because no 
corrections for the coupling between the nuclear and the 
electronic motion is introduced in this simplest form of 
the molecular Schrodinger equation. For the PEC V{R) 
the same analytical representation is used as e.g. in our 
earlier work on Ca2 [22j: 

It is divided into three parts: The central part of the 
potential (Ri < R < R a ) is given by 



V c (R)=T m + J2 
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The inner repulsive wall (R < Ri) is represented by 



(2) 



The long range part (R> R a ) is described by 
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V a (R) = U a 
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Here the aj are the fitting parameters. The radius R m 
forms together with T m the expansion point of the poten- 
tial and is typically chosen to be close to the equilibrium 
distance. Here it is taken as the radius of the minimum 
of a prior fitting result. The connection radii Ri and R a , 
as well as the parameters b and n are manually adjusted 
to get an acceptably small standard deviation using a low 
number of fit parameters etj . The application of the long 
range parameters slightly differs for the ground state and 
the excited state (see sec. [V]for details), the parameter 
C3 is needed only for the excited state, while the param- 
eters C6 to Cio are used for the ground state. 

Since the fit of this kind of potential description has a 
slow convergence, for the assignment of the data and to 
obtain RKR potentials as good starting potentials, a set 
of Dunham coefficients is derived first [23fl : 



88/88: v'=1,J'=1 63, v"= 
86/88: v -2, J - 24, v"= 
86/88: v'=2,J'=1 69, v"= 
86/88: v'=7, J - 99, v"= 
86/88: v'=7, J'=1 29, v"= 
88/88: v'=10,J'= 39, v"= 
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FIG. 1: An example spectrum averaged over 50 scans excited by a frequency of vl = 17397.62 cm -1 . In total, 14 progressions 
could be assigned which belong to the 2 1 £+ — X X E+ transition, the quantum numbers of only 6 of them are given on top of 
the figure including the isotope assignment. 



IV. ANALYSIS OF SPECTRA 
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The Yik are the Dunham parameters for the chosen 
reference isotopomer 88 Sr2 with reduced mass fiQ. T is 
the origin of the state. The reference of energy levels is 
defined by setting T for the ground state X*E+ to zero. 
The fa is the reduced mass of the isotopomer i, for which 
the energy ladder is considered. 



The main part of the analysis of data is done using a 
partially automated software, of which a prior version is 
already mentioned in [2J|. First attempts of assigning 
the spectra were done using the set of Dunham coeffi- 
cients reported in [4j together with a simple program, 
where the user has to select at least three lines from a 
possible progression, for which the program subsequently 
suggests assignments and automatically searches for ad- 
ditional lines within the spectrum. Using this function 
it became clear within short time, that the assumed ro- 
tational dependence of the term energies calculated with 
the Dunham set from [J| is incorrect even considering the 
relatively low accuracy of a few tenths of a wavenumbcr 
as given in [J]. In many cases the vibrational energy 



differences of odd rotational ground state levels of the 
reference isotopomer 88 Sr 2 were suggested as the ones 
fitting best, though these levels do not exist for homonu- 
clear isotopomers with nuclear spin zero like the even 
isotopes of Sr. In the diploma thesis by Schneider [25J . 
from which the paper |4| was derived, one can read that 
only the three strongest doublet progressions were used 
for the rovibrational assignment of the ground state and 
for the fitting of the Dunham-coefficient set which they 
used for the calculation of the RKR potential. Obviously 
their selected data set was too small to fix the assignment 
unambiguously, so our assignment procedure was started 
from the beginning. 



A. The initial assignment 

Initially, three progressions belonging to the same up- 
per v' and neighboring J' quantum numbers are taken, 
which were specifically produced by exciting at wave- 
lengths where collisionally induced rotational satellite 
lines appear. In this way the relative J and v assignment 
was fixed and only an offset in J of the whole set had to 
be varied to find the correct absolute rotational assign- 
ment. This was first done in a combined fit of Dunham 
coefficients for both electronic states, which was possible, 
because the rovibrational levels of the upper 2 1 £+ state 
do not show local perturbations in the small energy in- 
terval studied here. For the standard deviation a, which 
was minimized, we use the conventional definition: 
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Here N is the number of observations, M the number 
of fit parameters, and p n the uncertainty for the nth term 
energy or transition frequency. 

Among the first three measured progressions there was 
one which was already shown in [J] with the assignment 
(v' = 1, J' = 37). Varying the J offset of these three 
progressions yielded with almost the same quality three 
assignments, namely J' = 37, 39, or 41. After adding a 
fourth such progression to the data set, only an assign- 
ment to J'=39 or 41 remained possible, and after adding 
a fifth only for J' = 41 the standard deviation a was be- 
low 1.0 for an estimated uncertainty of p — 0.005 cm -1 
for each line. This rotational assignment then became 
unambiguous after adding more and more excitations to 
the data set. The data set was extended to higher J val- 
ues and to different upper vibrational levels by stepwise 
adding progressions. 

The vibrational assignment was taken from [J] and 
could later be confirmed by the observations of progres- 
sions belonging to the isotopomers 86 Sr 88 Sr and 87 Sr 88 Sr 
which were consistently described with the Dunham set 
by mass scaling (see eq. EJ. Also the results of the po- 
tential fits, which are described in Sec. IVl confirm this 



vibrational assignment. The v' assignment of the excited 
state is also confirmed by the intensity envelopes of the 
progressions e.g. originating from the lowest vibrational 
level of the 2 1 S+ state, i.e. v' = 0, the envelope shows a 
single maximum. Going to higher v' one sees as simple 
image of the wave function an envelope with the appro- 
priate number of maxima, v' + 1 . 



B. The automated assignment of the main body of 

data 



Besides the description of the ground state also an 
analysis of the excited state is developed. But each laser 
excited progression only gives information about one sin- 
gle rovibrational level of the excited state, so the informa- 
tion contained in the huge amount of weak lines belonging 
to overlapping excited progressions is very valuable, too. 
Since the amount of collisionally induced rotational satel- 
lite lines, which are intensively used, e.g. in \2A, [2&1 [27| . 
for the description of the excited states, is very low in this 
case, caused by the very weak overall intensity of the Sr 2 
fluorescence spectra, the analysis of these overlapping ex- 
citations is of special importance. Doing such an analysis 
by hand or even applying the half automated method as 
for the initial assignment of the strongest lines would 
have been very error prone and extremely time consum- 
ing, since it is not straightforward to decide which lines 
belong to the same progression. Even for a lot of com- 
binations of lines, which turn out later not to belong to 
the same progression, there is often more than one choice 
where they would fit nearly equally well into the ground 
state term energy set. 

To do the assignment with high reliability and with 
many consistency checks within limited time a spe- 
cial program code was written. The algorithm basi- 
cally checks for every unassigned line, starting with the 
strongest one, every possibility of assignment. First it is 
checked if the possible upper level could be excited by the 
actual laser setting or alternatively if it could be a colli- 
sionally induced satellite of an already assigned stronger 
progression. This information is collected for the present 
recording and the correctness of the assumed assignment 
is mainly judged by the number of lines found inside the 
spectrum which possibly belong to a common progres- 
sion, and in the case of a P/R doublet progression by 
the number of complete doublets. An important crite- 
rion is the maximal length of a progression, i.e. there is 
no assignment found which gives more Q lines or more 
complete P/R doublets in the recording. An additional 
acceptance criterion is the positive answer to the question 
if the upper level of the selected progression also fits into 
the simultaneously generated description of the excited 
state 2 1 S+. 

In the final assignment run of an iterative series of as- 
signments, of Dunham fits, and of potential fits (see sec- 
tion HIIli it was possible to identify around 10300 lines 
from the 37 recorded spectra within a computation time 
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FIG. 2: Overview of the observed energy levels of the ground 
state X 1 ^. 



FIG. 3: Observed energy levels of the excited state 2 S^j 



of about 8 hours. The resulting data confirm the applied 
model, because it describes them by the potentials and 
sets of Dunham coefficients with standard deviations a 
(see eq. ^ smaller than one. The uncertainties p are 
mostly four times smaller than the deviation limit set 
for the predicted value during the assignment, 0.02 cm -1 
for the ground state and 0.04 cm -1 for the excited state. 
Also very weak progressions of lines with signal-to-noise 
ratios down to 2.5 could be reliably assigned due to the 
many consistency checks. Only few single strong lines 
and some weak lines with S/N below 7 and very few 
medium strong lines with S/N below 30 remain unas- 
signed. The few unassigned strong lines form a doublet 
together with the laser line and are clearly identifiable as 
molecular lines because of rotational satellite structures, 
but do not fit into the Sr2 ground state potential as a P/R 
doublet. These lines possibly belong to another molecule 
as for example SrH or SrO. Few remaining unassigned 
lines could belong to different excited states, e.g. the 
lower 1 1 E+ state or the lowest two 1 H U states obtained 
by theoretical calculations [g, [[(J, LLL|], but which have, 
however, not been investigated so far. The main part of 
unassigned lines probably belongs to progressions which 
are too weak to match the requested assignment crite- 
rion of at least two clearly visible P/R doublets for each 
progression. 



C. Resulting data fields 

The data set obtained consists of about 4673 rovibra- 
tional energy levels for the ground state and 260 for the 
excited state 2 1 £+. For the ground state 3163 levels 
belong to the isotopomer 88 Sr 2 , 1059 to the isotopomer 
86 Sr 88 Sr and 451 to 87 Sr 88 Sr. The data fields are shown 
in figure[5]for the ground state and figure[3]for the excited 



state. 

Within a single progression the relative uncertainties 
of the transitions are estimated to be mainly about 0.005 
cm -1 from the applied resolution (0.05 cm -1 ) and the 
obtained S/N. The Doppler width is 0.03 cm -1 for a tem- 
perature of 950 °C and a frequency of 17250 cm -1 . Since 
overlapping excitations are evaluated, which contribute 
the by far largest amount of transitions to the data set, 
the excitation is not always on resonance and so lines of 
molecules with non zero velocity components parallel to 
the direction of observation are very likely. Thus Doppler 
shifts cannot be excluded. However, for the ground state 
only the differences of possible Doppler shifts of the indi- 
vidual lines within a progression contribute to the uncer- 
tainty. The depth of the ground state potential is below 
1100 cm" 1 , so the shift here is well below 0.005 cm -1 . 
The relative accuracy of the Fourier-transform spectrom- 
eter is 0.001 cm -1 . The errors of individual lines, which 
overlap with other transitions, are enlarged to 0.01 cm -1 . 

The term energies of the excited 2*E+ state are derived 
from the term energies of the ground state calculated 
from the X state potential (see section [Vj by adding the 
measured transition frequencies, and averaged over all 
observations of the same level. The uncertainties for the 
resulting energies are estimated to be 0.01 cm -1 . Since 
here the absolute transition frequencies are needed for 
determining the excited state levels the Doppler effect is 
the main source of error. Because of the overall weak- 
ness of the observed lines the probability of excitation 
in the wings of the Doppler-profile is rather low, so the 
estimation of an uncertainty of 0.01 cm -1 should be a 
conservative limit. 

The obtained standard deviations a (eq. O of the po- 
tential and the Dunham fits show that the given errors 
are rather overestimated than underestimated. 



TABLE I: Dunham coefficients for the X E+ ground state. 
All values are in cm -1 and for the reference isotopomer Sr2. 
For the range of quantum numbers see text. 
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V. RESULTS 
A. The X X S+ ground state 

The final set of Dunham coefficients with 35 freely vari- 
able parameters for the ground state is constructed by 
fitting the calculated level differences to differences of 
transition frequencies within the individual progressions. 
The result can be found in table Q] and describes those 
rovibrational levels shown in fig. [2] which are on the left 
side of a line defined by the points (v"=47, J"=79) and 
(v"=17, J"=219). The standard deviation is a = 0.77. 
For the description of the asymptotic energy levels by 
Dunham coefficients an inappropriately large amount of 
coefficients would be necessary and even a convergence is 
not guaranteed. The Dunham parameters should be used 
with caution, a physical meaning can only be attributed 
to the lowest ones like Y\o, loi and lo2- 

For the fit of the X state potential curve (see section 
Mil for the description of the analytical model) the ener- 
gies of the originating excited levels of the progressions 
are handled as free parameters. In this fit the energies 
are calculated with reference to the ground state asymp- 
tote. Therefore, the potential parameter U^, is set to 
and the parameter T m is adjusted to get a continuous 
connection of the central potential part (Ri < R < R a ) 
to the long range area (R> R a )- The parameter C3 is 
for the ground state while the parameters Cq and Cg are 
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FIG. 4: The residuals of the X 1 E^" potential fit as a function 
of the transition frequencies. The range S marks the uncer- 
tainty interval estimated for most lines. 



taken from theory [28| . To improve the description of the 
data the parameter Ciq (also theoretical values exist) was 
varied. The parameter A is adjusted to get a continuous 
connection to the inner potential wall (i? < Ri, formula 
([3])), while the parameter B is varied in the potential fit. 

Figure 0] shows the residuals of the potential fit as a 
function of the energy, while the resulting potential coef- 
ficients (only 22 free parameters compared to 35 Dunham 
coefficients) are given in table [JT] The weighted standard 
deviation of the fit is a — 0.81. The potential was ap- 
plied for the integration of the Schrodinger equation in 
the interval from 3.3 A to 100 A. 

To improve the description of the long ra nge par t we 
also used information from trap experiments [2ll |29( de- 
rived by photoassociation. In these experiments the pho- 
toassociation laser was tuned to excite high vibrational 
levels of the 2 1 S+ state. Such resonances lead to trap 
losses the magnitude of which is directly related to the 
Franck-Condon density (the square of the overlap integral 
of the scattering wave function of the ground state and 
the vibrational wave function of the 2 1 S+ state) of the 
excited transitions. Because the vibrational wave func- 
tions of bound states have significant amplitudes mainly 
close to the outer turning points, the intensity envelope of 
the photoassociation resonances gives information about 
the amplitude of the scattering wave function and thus 
about the node positions in this region. Here the posi- 
tion of the second to last node for the molecule 86 Sr2 was 
obtained by [2l[ to be at 62.6 ao and the position of the 
last node for the molecule 88 Sr2 slightly below 40 A by 

To include this information in our potential determi- 
nation, we first adjusted the C10 coefficient to fit the po- 
sition of the second to last node of the isotopomer 86 Sr2 
to the value above and calculated the term energy for the 
last vibrational level (J"=0, v"=62) with this potential. 



TABLE II: Potential coefficients for the X 1 E^" ground state, 
referenced to the asymptote. 
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cm -1 


b 




-0.57 


r^m 


4.67169686 A 


T 


-1081.8163 


cm 


R* 


3.98 A 


ii 




6 


A 


-2.5412398xl0 3 


cm -1 


B 


9.88561567xl0 6 cm" 1 A 6 


Ho 


11.0 A 


C 6 [28] 


1.4955xl0 7 cm -1 A 6 


C 8 [28] 


5.1175x10 s cm" 1 A 8 


Cio 


2.495 xlO 10 cm 


-i A io 


Uoo 


0.0 


cm" 


derived constants: 






D e 


1081.82 


cm -1 


Rc 


4.67174 A 



This term value was included in a new overall fit as an 
observed quantity with an error of 3 x 10" 5 cm" 1 . Fi- 
nally, it was checked that this potential also reproduced 
quite well the position of the last node of the scattering 
wave function of 88 Sr2 reported by Yasuda [291 ]. 

Figure [5] shows the final potential, which is valid in 
the region from 4 A to 11 A for our spectroscopic data, 
compared to the RKR potential [J| and the theoretical 
calculations [g]. Here it is clearly visible, that the devi- 
ations of the two experimental potentials caused by the 
change in the rotational assignment are mostly bigger 
than the differences to the theoretical calculations. 




R[A] 



FIG. 5: The X 1 E+ ground state potential (solid line) com- 
pared to the RKR-potential from 4] (dotted) and the ab initio 
potential from [Q] (dashed dotted). 



TABLE III: Dunham coefficients for the excited state 2 Ej, 
for the isotopomer 88 Sr2. All values in cm" 1 . T is referenced 
to the origin of the Dunham series for the ground state. 



Ufe-> 





1 


2 








0.0219691 


-5.962xl0" 9 


i 


84.207 


-6.693xl0" 5 


-6.17xl0 -11 


2 


-0.2639 


-6.38xl0" 7 




3 


-0.001422 






T 


17358.7262 







B. The excited 2 X E+ state 

Table IIIII shows the set of Dunham coefficients with 
9 freely varied parameters derived in a combined fit for 
both electronic states from 260 levels of the low part of 
the 2 1 E+ state. It reproduces all observed energy levels 
with v' up to 12 and J from 9 to 221 with a standard 
deviation of a = 0.80. 

For the final potential description of the excited state 
2 1 E+ term energies are derived, which are referenced to 
the potential minimum of the ground state. As already 
mentioned in section IIVCI these term energies are cal- 
culated by adding the measured transition frequencies 
to the term energies calculated from the ground state 
potential. The model description differs slightly from 
that of the ground state. Here the parameter T m is in- 
cluded in the fit while the parameter U&3 is calculated by 
adding the ground state dissociation energy of 1081.82 
cm" 1 to the atomic transition frequency for 1 So — ► Pi 
of 21698.482 cm" 1 obtained from [30J]. The long range 
coefficient C3 is taken from [29| , while the Cq and Cg co- 
efficients are adjusted to get a continuously differentiable 
connection of the long range part (R > R a ) to the cen- 



TABLE IV: Potential coefficients for the excited state 2 J E+ 
with reference to the minimum of the X state. 



ai 


1.619x10° cm -1 


a 2 


2.0176667xl0 4 cm" 1 


a 3 


1. 049518x10* cm -1 


EI4 


-1.47839xl0 4 cm" 1 


as 


-4.65051 xlO 4 cm -1 


as 


-5.36901 xlO 4 cm -1 


b 


-0.50 


Rm 


4.1783479 A 


-L m 


17358.7389 cm" 1 


R. 


3.775 A 


11 


6 


A 


1.4855497xl0 4 cm -1 


B 


1.01060899xl0 7 cm -1 A 6 


Ra 


5.0 A 


C 3 [29] 


5.9712xl0 5 cm _1 A 3 


C 6 a 


-5.1541xl0 7 cm" 1 A 6 


C 8 ° 


8.9977xl0 8 cm -1 A 8 


Uoc 


22780.30 cm -1 


derived constants: 




T e 


17358.739 cm" 1 


Re 


4.17828 A 



"The coefficients Cg and Cg are only parameters for connection 
to the atomic asymptote 1 So+ 1 Pi. 

















ae Sr 




















0.015- 










■ 


86 Sr 88 Sr 










.o" . 

o 


. ■ . 


o 


87 Sr 88 Sr 




















0.010- 




















■ 




" ■ " T) . "" -P 


. ■ 




" 






















'f 


0.005- 


■ 




p: % . 


■ ■' ■ 


■ 






c^ 






. 


■%. o ■ '. ■ ■ .■ 










_U) 


0.000- 


¥ 




■ . ■ : ">■■■■■ 


■ o 


m m 




8 


"O 

tn 

0) 


-0.005 - 


% 




■ v>v.-.':-° 

■rita , J ■ ■ ■ ■ 

* 1 ■ "." , 


■ 


■ 






cn 


-0.010- 


■ "■ 




■'.■. :■*■■*■.. 

. * ■ 


■■ 




























. 




■ .- ■■■ 










-0.015- 















17400 17600 17800 18000 18200 18400 18600 18800 
Term energy [cm 1 ] 



FIG. 6: The residuals of the 2 Ej potential fit as a function 
of the term energies. 



tral region (R4 < R < R a ). The parameters A and B are 
adjusted to get a continuously differentiable connection 
of the inner potential wall (R < Ri). 

The 2 1 £+ potential supports roughly 350 bound vi- 
brational levels, with a large uncertainty, ± 50 or more, 
because t he gap b etween our data and the photoassocia- 
tion data [29|, |3l| is about 4000 cm" 1 . 




3.8 4.0 



4.2 4.4 

R[A] 



FIG. 7: The derived potential (solid line) compared to the 
RKR potential from Gerber [|[ (dotted) for the state 2 1 E+ 
and the ab initio potentials from Boutassetta [a, [33] (dashed 
dotted) and Czuchaj [lOi. 1331] (dashed). 



The residuals of the potential fit, which has a weighted 
standard deviation of a = 0.68 (defined in eq. ©), can 
be found in figure [51 The potential coefficients are listed 
in table HVl The integration interval used for solving the 
radial Schrodinger equation extends again from 3.3 to 
100 A. 

Figure gives a comparison of the resulting potential 
energy curve (PEC) to the experimental potential from 
4] and to the theoretical calculations [6[ and [lCJ]. Here 
the large deviations between the two experimental po- 
tentials are again obvious, due to the different rotational 
assignments of the excited state, which are directly re- 
lated to the ground state. 



VI. CONCLUSION AND OUTLOOK 

In this article two electronic states of the Sr2 dimcr 
were studied, the ground state X X S+ and the minimum 
region of the excited state 2 1 E+. During the analysis 
the rotational assignment had to be shifted by four units 
upwards compared to the previously most precise investi- 
gation of this molecule [4| . This also implies a significant 
shift of the potentials towards larger inter atomic dis- 
tances. The observed 3163 ground state levels for the 
isotopomer 88 Sr2 are more than half of the existing rovi- 
brational bound levels of this isotopomer (around 5400). 
The observed ground state levels of the isotopomers 
86 Sr 88 Sr (1059) and 87 Sr 88 Sr (451) further confirm the 
new rotational and the earlier vibrational assignment. 

We improved the asymptotic range of the ground state 
potential by the use of photoassociation data [2l|, [29|, [31( 
and calculated the scattering lengths for all combinations 
of natural Strontium isotopes. However, there remains 



TABLE V: The scattering lengths for the different isotopic 
combinations calculated using the long range coefficients from 
[341 ] and from [28|], all values in atomic units (ao = 0.5292 A). 

isotopomer fit using [34] fit using [281 

127 
36 
-45 
-1900 a 
1430 
171 
101 
99 
58 
4.5 



84 Sr+ 84 Sr 


121 


84 Sr+ 86 Sr 


30 


84 Sr+ 87 Sr 


-64 


84 Sr+ 88 Sr 


1170 


86 Sr+ 86 Sr 


677 


86 Sr+ 87 Sr 


160 


86 Sr+ 88 Sr 


97 


87 Sr+ 87 Sr 


95 


87 Sr+ 88 Sr 


54 


88 Sr+ 88 Sr 


-4.8 



a For the isotopic combination 84 Sr+ 88 Sr the number of vibra- 
tional levels decreases from 63 to 62 when the long range coefficient 
set from 28] is used instead of that from 34] , thus only a large mag- 
nitude is predicted here. Also the sign in the case of 88 Sr+ 88 Sr is 
not determined, while in this case the scattering length is close to 
zero. 



still a significant uncertainty which stems from the un- 
certainty in the long range coefficients, especially from 
that of the C 6 coefficient, as discussed in [2l|. Regard- 
ing this coefficient five different theoretical calculations 
[23, H, [33, ys [33] were published during the recent years, 
with values ranging from 3103 a.u. [28| to 3249 a.u. [34| . 
Since just in these two publications with results deviat- 
ing strongest the authors have used the most advanced 
methods and are the only ones who give also Cs and C10 
coefficients, we performed two different potential fits ap- 
plyin g bo th sets of Q. The first one using the coefficients 
from [28| is shown in table [TT1 

The dissociation energy D P differs for both fits: Us- 
ing the coefficients from [28f we obtained D P = 1081.82 



cm , and with those from 



D, 



1081.52 cm- 



In 



fits to reproduce the scattering node position from [21( 
for the potential using the coefficients from [3J] the co- 
efficient C10 had to be reduced from 4.25 x 10 7 a.u. to 
1.52 x 10 7 a.u. while for the potential using the coef- 
ficients [28[ it had to be increased from 4.22 x 10 7 a.u. 
to 6.60 x 10 7 a.u. We have no argument which solution 
should be preferred, thus we need more measurements. 
The number of vibrational levels, 62 for the three lighter 
isotopomers 84 Sr2, 86 Sr 84 Sr and 87 Sr 84 Sr and 63 for the 
six heavier isotopomers, is the same in both cases, only 
for 84 Sr 88 Sr it changes from 62 for the long range co- 
efficients from 34] to 63 for the coefficients from [28J J . 
respectively. 

The resulting scattering lengths are compared in ta- 
ble [Vj Here the two results for the isotopic combination 
86 Sr+ 86 Sr are clearly located inside the range given in 
[211 ] (610 ao to 2300 ao), as expected because of the use 
of their node position of the scattering wave function, 



but for the combination 88 Sr+ 88 Sr the value of -4.8 ao 
calculated using the long range coefficients from [34j is 
outside the range of -1 ao to 13 ao reported in [2l| . Mick- 
elson calculated |21| the scattering length for the isotopic 
combination 88 Sr+ 88 Sr by mass scaling from the result 
for the combination 86 Sr+ 86 Sr using either the ab initio 
potential [6| or the RKR potential Q, which both sup- 
port a lower number of vibrational levels compared to 
our result. 

Another uncertainty for the prediction of the scatter- 
ing lengths lies in the accuracy of the Born-Oppenheimer 
approximation for this molecule, though possible correc- 
tions necessary for mass scaling should be small because 
of the large mass of the atoms. 

The derived potentials present an excellent starting 
point for further spectroscopic investigations of this 
molecule, concentrating on a more complete investiga- 
tion of the long range part of the ground state. This 
will help to reduce the uncertainties in the calculation 
of the scattering lengths. The calculated Franck-Condon 
factors for the transition 2 1 S+ -- X 1 !]"!" direct us that 
the highest vibrational levels for 88 Sr2 could be observed 
through fluorescence by laser excitation of the v'=4 level 
of the state 2 1 E+ starting from a low vibrational level 
of the ground state as e.g. v" = 2. An investigation 
is presently prepared. If this should not give satisfy- 
ing results, a different technique will be applied, where 
the v'=4 level becomes excited directly from the desired 
asymptotic levels, which are significantly thermally pop- 
ulated due to the high temperatures and the flat ground 
state potential. Because of the high density of vibrational 
levels close to the ground state asymptote the number of 
overlapping bands will be very high, too. A monochro- 
mator could be used to filter the fluorescence in the range 
where only transitions from the v'=4 level to a certain 
low ground state level are expected, for example to the 
frequency of the transition v' = 4 — v" = 0. This way 
the selective observation of the desired transitions was 
already successfully applied for the Ca2 dimer [3a |. 

The transition path for reaching asymptotic levels from 
v"=0 or 2 via v'=4 can be inverted to produce ultra- 
cold molecules in v"=0 from photoassociation followed by 
spontaneous decay for populating the asymptotic level. 
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